Antiferromagnetism and d-wave superconductivity in cuprates: a cluster DMFT 

study. 
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We present a new approach to investigate the coexistence of antiferromagnetism and d-wave su- 
perconductivity in the two dimensional extended Hubbard model within a numerically exact cluster 
dynamical mean-field approximation. Self-consistent solutions with two non-zero order parameters 
exists in the wide range of doping level and temperatures. A linearized equation for energy spectrum 
near the Fermi level have been solved. The resulting d-wave gap has the correct magnitude and 
k-dependence but some distortion compare to the pure A x i_ y i superconducting order parameter 
due to the presence of underlying antiferromagnetic ordering. 
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A microscopic theory of High-temperature supercon- 
ducting cuprates (HTSC) is still far from the final un- 
derstanding One of the most important recent 
experimental achievements was the discovery of pseudo- 
gap (PG) phenomenon above superconducting transition 
temperatures (T c ) 0] and existence of a sharp 41-meV 
resonance below T c related with some collective antifer- 
romagnetic excitations Q . Recent neutron scattering ex- 
periments Q provide a new insight for the interesting 
problem on the origin of a condensation energy. The co- 
existence of an antiferromagnetism (AFM) and d-wave 
superconductivity (d-SC) in cuprate could be a natural 
way of discussing such different HTSC phenomena. This 
require a quantitative electronic structure theory includ- 
ing two different type of the order parameters: AFM and 
d-SC. Within such approach one can in principle ana- 
lyzed the phase diagram of HTSC cuprate and interplay 
between antiferromagnetism and d-wave superconductiv- 

ity §■ 

A minimal theoretical tool for cuprates consists of the 
two-dimensional Hubbard model Q|. The importance of 
including the realistic tight-binding spectrum obtained 
from the LDA-band structure analyses M was realized 
during the last years. Unfortunately, a most accurate 
Quantum Monte-Carlo (QMC) simulation of hole-doped 
2d-Hubbard model has a difficulty to describe an inter- 
esting part of the HTSC phase diagram near 15% doping 
at the low temperature due to so-called sign-problem || . 
The perturbation theory of d-SC |[(| ignore the vertex 
corrections in the strong correlation case of HTSC. Great 
progress in the theory of the interacting fermions results 
from the developing of the dynamical mean-field theory 
p|,^2[ . While the antiferromagnetic phase is easy to in- 
corporate in the single-impurity DMFT-approach p2[ , 
the d-wave superconductivity requires a cluster general- 
ization of DMFT. Different cluster-DMFT scheme have 
been proposed [jl2],[l3| and the recent application to the 
problem of the pseudogap in HTSC [H have shown the 
efficiency of the cluster DMFT approach. The investiga- 
tion of paramagnetic phase of two dimensional Hubbard 



model could be simplified using a translational symmetry 
p3| , while the problem of a coexistence of AFM and dSC 
demands a broken-symmetry cluster calculation. This 
is equivalent to multi-orbital DMFT approach EE] and 
could be solved with QMC method @. 

In this Letter we investigate the problem of anti- 
ferromagnetism and d-wave superconductivity in two- 
dimensional Hubbard model within a cluster DMFT 
scheme. 

The minimal cluster which allow us to study both AFM 
and d-SC order parameters on the equal footing consists 
of 2x2 system in the effective DMFT-medium. We start 
with the extended Hubbard model on the square lattice: 
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where ty is an effective hopping and Ui local Coulomb 
interactions. We chose nearest-neighbor hopping t = 
0.25 eV and the next nearest hopping t'/t = —0.15 
for the model of La 2 - a; Sr ;c Cu04 ||. The total band 
width is W=8t and the Coulomb interactions set to be 
U/W=0.6. Let us introduce the "super-site" as an 2x2 
square (Fig.|l|). The numeration of the atoms in the 
super-site is also shown in the Fig.[l]. It is useful to in- 
troduce the superspinor = {cf a } where a = 1, 2, 3, 4. 
(the spin indices are not shown) . Taking into account the 
spin degrees of freedom, this is 8-component superspinor 
creation operator. Then the crystal Green function for 
the Hubbard model can be rewritten as 

G (k,iw) = [iu> + n — h (k,iu;)] 1 

where h (k,zw) is the effective hopping supermatrix and 
H is the chemical potential. For simplicity we will write 
all the formulas in the nearest-neighbor approximations: 
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2x2 matrix in the spin space. Within cluster-DMFT ap 
proach we introduce intraatomic self-energy So and in- 
teratomic self-energies £, 
in the sense of the super-site: 



E y and the both are intra-site 
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The effective Hamiltonian defined through the transla- 
tionaly invariant self-energy corresponds to t x ~ t + E x , 
t y = t + E y renormalized energy dependent hopping. 
The functions Eo (ito) , E^ (ioj) , E^fiw) are found self- 
consistently within DNFT scheme |12[ and for the d-wave 
superconduction state E K ^ Ej,. It is straightforward to 
generalize this scheme for next-nearest neighbor (or more 
extended) hopping as well as the long-range Green func- 
tion and self-energy. In this case we can renormalized also 
the second-nearest hopping: t xy = t' + T, xy for the 2x2 
cluster, where E xy (or E 02 ) is the non-local self-energy 
in xy direction. 

According to the prescription of DMFT scheme Jl2| , 
we can write the matrix equation for the so-called bath 
Green function matrix Q which account for a double 
counting correction for the self-energy: 

g^iiiu) = G" 1 (iui) + E {iui) 

where the local cluster Green function matrix is equal to 
G a p (iui) = G a p (k,iw) , and summation is run over 

k 

the Brillouin zone of the square lattice. Note that in 
theEq.(g) we use translationaly invariant self-energy ob- 
tained from the cluster DMFT similar to the dynamical 
cluster approximation p3| . The present "matrix" form of 
a cluster DMFT with the self-energy which is not periodic 
inside the cluster allow us to study a multicomponent or- 
dered state. 

In this case we have the standard DMFT problem with 
four "orbital" states per super-site. It is solved by multi- 



orbital QMC technique [116 



For the problem of a coex- 



istence of magnetic ordering and superconductivity one 
can use the generalized Nambu technique |l7]] . We intro- 
duce the superspinor 



and the anomalous averages describing the (collinear) 
antifcrromagnctism (cf^Cji) and the superconductivity 
= (c a c iT ). 

We use the generalization of the Hirsch-Fye QMC- 
algorithm jOl for superconducting problem jl9|. In the 
4-spinor case a discrete Hubbard-Stratonovich transfor- 
mation has the following form: 
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where Xi = iarccosh[exp(i AtL^ 
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FIG. 1. (a) Schematic representation of antiferromagnetic 
d-wave 2x2 periodically repeated cluster; (b) generic phase 
diagram of HTSC materials; (c) The calculated values of two 
order parameters (local magnetic moment and d-SC equal 
time Green function F 01 ) for different hole doping (x) at the 
inverse temperature /3 = 60 eV _1 (T = 190K). 

We will take into account only the singlet pairing and 
in the case of d-wave there are following nonzero ele- 
ments of A matrix: A = A12 = — A23 = A34 = — A41. 
One can chose Ay to be real and therefore symmetric: 
Aij = Aji. Separating normal and anomalous parts of 
the Green function we have 



G(k,r,r') 



G(k,T,r') F(k,T,r>) 
F+(k,r,r') -G(-k,r',r) 



where G (k,r, r') = - (T r G k (r) C+ (r')) , F (k,r, r>) = 
— (T T Ck (r) C_k (t')) are the matrices in spin and "or- 
bital" space. It is convenient to expand the anomalous 
Green function in Pauli matrices F — (F° + Fa) ia v and 
use the symmetry properties [Q : 

F°(k,r,r') =F°(-k,r',r) 
F(k, r,r r ) =-F(-k,r',r) 

Therefore in the collinear antiferromagnetic case with d- 
wave superconductivity the 4x4 spinor formalism is re- 
duced to 2 x 2 one with the following spin-matrix form 
of the local Green function for our super-site: 



G(t,t') 



G T (r, 



F(t,t') 
-Gi (t',t) 



and the QMC formalism for the antiferromagnetic su- 
perconducting state is equivalent to the previous non- 
magnetic one |l9[ ]. Using the discretization of [0,(3] in- 
terval with L-time slices (At = (3/L and (3 — l/T the 
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inverse temperature) G a and F Greens functions are ma- 
trices of the dimension 2NL, where ./V is the number of 
atoms in the cluster. After a Fourier transform to the 
Matsubara frequencies we have: 



G (iu) = 



G T (iu) F (iu>) 



F(iu>) -G\ (iu) 
In superconducting states JlJ] the self-energy defined as: 

and the crystal Green functions is equal to: 



G" 1 (k,^)= ( ^ + M -M k ^) 
v ' V s (k,zw) 



s (k,iw) 
iuj — fj, + h* (k,io;) 



where s (k,io>) translationaly invariant anomalous part of 
self-energy S(iuj) similar to Eq. ([!]). 
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FIG. 2. Imaginary time normal (G) and superconduct- 
ing (F) Green functions for the 2x2 cluster DMFT solution 
with second-nearest neighbor hopping and inverse tempera- 
ture 13 = 50 (eV- 1 ) (T=230 K). 

The two-component order parameters state which con- 
sist of Neel antiferromagnetic state and d-wave super- 
conducting (Fig. la) lowered the symmetry of effective 
cluster-DMFT problem and did not allowed to use single- 
atom translational symmetry. The self-consistent DMFT 
8x8 matrix cluster problem with AFM and d-SC gen- 
eral order parameters have been solved within QMC 
scheme with L=64 time slices. The resulting two or- 
der parameters for /3 = 60 (or T=190 K) and t' = 
was presented in Fig.lc together with the generic HTSC 
phase diagram (Fig. lb) as function of the hole dop- 
ing. In this case the ordered magnetic moment is di- 
rectly related with imaginary-time Green function G ct (t): 
M = Gf (0) - (0) and for d-SC order parameter 
we chose a positive value of superconducting imaginary 



time Green function F 01 (0). It is interesting that the 
AFM cluster-DMFT solution exist for a much higher dop- 
ing concentration then experimental AFM ordered state 
and describe a dynamical mean-field version of AFM-spin 
fluctuations related to pseudogap phenomena (PG-region 
on the fig. lb). The maximum of d-SC order parameter 
corresponds to the doping level of about 15% and agree 
well with the generic HTSC phase diagram. The d-SC or- 
der parameter is zero close to the undoped region (x=0), 
due to presence of a large AFM-gap. When the AFM- 
gap is closed (x~ 5%) the d-SC states developed and after 
x~ 20% it decreasing again since AFM spin-fluctuations 
around (7r,7r) point disappear The precise charac- 
teristic of the phase diagram including the interactions 
between AFM and d-SC order parameters demands an 
extensive cluster-DMFT calculations for different tem- 
peratures and doping. 

We would like to note that the existing of d-SC cluster- 
DMFT solution for such high temperature does not nec- 
essary means that the superconducting transition tem- 
perature is even larger then 190K in our model. A crude 
estimation shows that d-SC solution disappear about 300 
K for x=0.15 and AFM solutions for x=0 become unsta- 
ble at the temperature just above 1000 K. This could be 
the sign of a "local" AFM solutions or local d-wave so- 
lutions, like a local moment in magnetic systems. We 
plane to estimate a proper superconductiong transition 
temperature in the future publication. It is important, 
that we find no serious sign-problem for all QMC calcu- 
lations with various doping level, probably due to "sta- 
bilized" antiferromagnetic dynamical mean fields acting 
on the neighbor atoms to our 2x2 cluster. 




FIG. 3. The d-wave gap function at the Fermi surface for 
x=0.15, t' = and /3 = 50 (eV -1 ). 
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The general role of next-nearest hopping is the low- 
ering of the van-Hove singularity Q which increase the 
density of state at the Fermi level for the doped case and 
favored d-SC solution. There is also change in the spin- 
fluctuation spectrum related with broadening of AFM- 
peak near (jr, tt) point due to formation of so-called ex- 
tended van-Hove singularities with increasing of the t'. 
On the Fig.| we show one of the AFM-dSC solution 
with the next nearest-neighbor hopping for the 10% dop- 
ing level and (3 = 50 (ev~l). Resulting local magnetic 
moment is M=0.28 [i B and F(0)=0.036.Onc could see 
that the superconduction order parameter is really of 
the d x 2_ y 2 symmetry since diagonal elements (F 00 ) as 
well as the nearest-neighbors elements (F 00 ) are all equal 
to zero and the only nearest neighbor superconducting 
Green function (F 01 ) is non-zero and moreover change 
the sign for F^ and F y components. The normal local 
Green function (G 00 ) (ploted for the spin-up atom in 
the Fig.||) as well as (G 02 ) is spin-split, while nearest- 
neighbor Green function (G 01 ) has no spin-splitting due 
to AFM spin symmetry (see Fig.|l|). 

In order to find the superconducting energy gap we 
solved the linearized equation for energy spectrum, as- 
suming that the characteristic energy scale of E CT (iio) 
and S (iu>) is larger than the SC-gap A (icu). In this case 
we could perform analytical energy continuations and the 
generalized equation for the energy spectrum has the sim- 
ple form: 

det(i? - EO) = 

where H = t (k) + S(0) - fi , O = 1 - S'(0) and 

S(0) = Jjf E(r)dr , E'(0) = $ rS(r)dr. Note that 
E(0) and S'(0) in this expression should be also trans- 
lationaly invariant similar to Eq. ([!]). We solve the lin- 
earized equation for energy spectrum (for t' = and 
(3 = 50) and obtained the superconducting energy gap 
on the Fermi surface (Fig.|J). The topology of the Fermi 
surface was defined as the zero-energy contour for the 
energy spectrum with all F's Green functions equal to 
zero. It is clear that symmetry of the d-wave state is 
not pure d x 2_ y 2 due to underlying AFM ordered states 
which lowered the symmetry of d-SC. This also means 
that d-SC order could lower the symmetry of the AFM 
Neel state and more general non-collinear magnetic state 
need to be investigated. Nevertheless the gap function 
has a maximum near the (tt, 0) and (0, ir) points and is 
almost zero near the (tt/2, it/2) point and the magnitude 
of the maximum superconducting gap is of the order of 
15 meV in a good agreement with experimental estima- 
teess Q and much lower then the AFM gap for undoped 
case. 



In conclusion, we present the evidence for coexisting 
of antiferromagnetism with the d-wave superconducting 
states within the cluster DMFT two dimensional ex- 
tended Hubbard model. 

We thank Gabi Kotliar, Antoine Georges and Yurii 
Izyumov for helpful discussion. 
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